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ABSTRACT 

We present ultra-deep J and Ks imaging observations covering a 30' x 30' area of the Extended 
Chandra Deep Field-South (ECDFS) carried out by our Taiwan ECDFS Near-Infrared Survey (TE- 
NIS). The median 5a limiting magnitudes for all detected objects in the ECDFS reach 24.5 and 23.9 
mag (AB) for J and Ks, respectively. In the inner 400 arcmin 2 region where the sensitivity is more 
uniform, objects as faint as 25.6 and 25.0 mag are detected at 5a. So this is by far the deepest J 
and Ks datasets available for the ECDFS. To combine the TENIS with the Spitzer IRAC data for 
obtaining better spectral energy distributions of high-redshift objects, we developed a novel decon- 
volution technique (IRACLEAN) to accurately estimate the IRAC fluxes. IRACLEAN can minimize 
the effect of blending in the IRAC images caused by the large point-spread functions and reduce 
the confusion noise. We applied IRACLEAN to the images from the Spitzer IRAC/MUSYC Public 
Legacy in the ECDFS survey (SIMPLE) and generated a J+Ks selected multi-wavelength catalog 
including the photometry of both the TENIS near-infrared and the SIMPLE IRAC data. We publicly 
release the data products derived from this work, including the J and Ks images and the J+Ks 
selected multiwavelength catalog. 

Subject headings: catalogs — cosmology: observations — galaxies: evolution — galaxies: formation 
— galaxies: high-redshift — infrared: galaxies 



1. INTRODUCTION 

Near-Infrared (NIR) imaging surveys provide several 
advantages over optical observations for studies of galax- 
ies. For nearby galaxies, stellar light in the NIR better 
traces the dominant stellar population by mass and is less 
affected by dust extinction. For distant galaxies, either 
the Balmer break is shifted to the NIR bands (ellipticals 
at z > 1.5) or the optical flux may be obscured by dust 
(dusty starburst galaxies or AGNs). These effects pro- 
duced red optical-to-NIR colors, and at their most severe 
extremes result in var ious red galaxy populations ( e.g., 



extremely red obj ects. lElston. Rieke. fc R ieke 1988; dis- 
tant r ed galaxies. iFranx et al.]2003l: s ee also. Yan et al 
2004 iWang. Cowie. fc Barger I 120121 and IGuo et al. 



2012ft . The most distant Lyman Break Galaxies (LBGs) 
known (z > 7) cannot be detected at all in most op- 
tical bands because of absorption by the inter-galactic 
medium, and can only be detected in the NIR. Apart 
from studies of galaxies, NIR observations also are impor- 
tant for studies of Galactic dwarf stars and starforming 
regions (both are very red because of cither low surface 
temperature or dust extinction). Hence NIR imaging 
surveys are very valuable for studying Galactic objects, 
stellar masses and the distribution of the dominant stel- 
lar population by mass in nearby galaxies, and for identi- 
fying as well as studying the star-forming properties and 
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the evolved population of distant galaxies. 

The M ultiwavelength Surve y by Yal e-Chile 
CMUS YC . IGawiser et al~l \xM iTavlor et all l200l 
iCardamone et al. 112010ft covers the 30' x 30' Extended 
Chandra Deep Field-South (ECDFS) with observations 
in optical and NIR. The depths of the NIR data from the 
MUSYC observations, however, are too shallow to study 
the most distant galaxies, (23.0 mag and 22.3 mag for J 
and K , respectively, at 5a for point sources) and are not 
comparable to those of the Spitzer IRAC observations 
performed by the Spitzer IRAC/ MUSYC Public Legac y 
in the ECDFS survey ( SIMPLE, [bamen et all 12011ft . 
According to the curr ent studies of the luminosity 
funct ion at z ~ 7 (e.g. lOuchi et al. 1 120091 lYan et alH 
2011), an L* galaxy at z = 7 would have apparent AB 

magnitudes of 26 27 in J and Ks, and about one 

galaxy with J ~ 25 at z ~ 7 is expected to be found 
in a field size similar to that of the ECDFS. To find 
LBGs at z > 7 for constrainin g the very bright e nd of 
the z > 7 luminosity function (jHsieh et al. 1 12012ft and 
to study the properties of dusty galaxies at z = 2-5, 
we initiated the Taiwan ECDFS Near-Infrared Survey 
(TENIS) in 2007. This survey comprises extremely deep 
J and Ks imaging observations in the EC DFS using the 
Wide-field InfraRed Camera (WIRCam, iPuget et al. I 
12004ft on the Canada-France- Hawaii Telescope (CFHT). 
The 5a limiting magnitudes for point sources achieve 
25.6 and 25.0 mag in J and Ks, respectively, which 
shows that the TENIS data are by far the deepest NIR 
data in the ECDFS. 

The SIMPLE project provides deep IRAC data at 
the wavelengths of 3.6, 4.5, 5.8, and 8.0/im in the 
ECDFS. In addition, in the central 10' x 15' area of the 
ECDFS, the Great Observatories Origins Deep Survey- 
South (GOODS-S, IGiavalisco et~aT1 l200l project pro- 
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vides ultra-deep IRAC data (M. Dickinson et al., 2012, in 
preparation). These IRAC data are very important for 
various researches especially for high-z studies. There 
are many continuum features of distant galaxies shifted 
to wavelengths beyond the K$ band. For example, for 
galaxies at z > 0.4, the rest-frame 1.6/Ltm bump caused 
by the H~ opacity minimum in the stellar ph otosphere 
([Simpson fe Eisenhardt 1 119991: iSawicki 1 12002ft and the 
Balmer break for galaxies at z > 4.5 are redshifted out 
of the Ks band. Properly detecting these features in 
IRAC bands can improve the quality of the photomet- 
ric redshift estimation. In addition, the Balmer break is 
also an important tracer of ag e and stellar mass o f high- 
redshift galaxies. Moreover, iHsieh et al. I (|2012ft show 
that IRAC data are essential for minimizing contamina- 
tion from Galactic cool stars when search ing for z > 7 
LBGs, and lWang. Cowie. fc Barger I ((2012ft show that K 
- IRAC color is able to pick up the most extremely dust- 
hidden galaxies at redshifts between 1.5 and 5. These 
studies show the importance of IRAC data for high-z 
studies. Measurements of IRAC fluxes for individual 
sources, however, can be easily contaminated by blended 
neighbors because of the large IRAC point-spread func- 
tions (PSFs). This issue becomes very serious in deep 
IRAC surveys (e.g., the GOODS IRAC survey) where 
the surface density of faint objects becomes very high. To 
improve the ECDFS IRAC flux measurements, we have 
developed a novel deconvolution technique to reduce the 
effect of object blending in the SIMPLE images. We 
then combined our J, Ks data with IRAC photometry 
to form a multiwavelength catalog. 

We present the survey fields of the TENIS, SIMPLE, 
and GOODS-S IRAC projects in Figure [TJ We describe 
the observations of TENIS in Section [2J the procedure 
of the TENIS data reduction in Section and the 
evaluation of the reduction quality in Section [U We 
describe the photometric catalog of TENIS J and 
K$ data in Section [5] and compare it with other NIR 
datasets in Section [6] In Section [71 we present our newly 
developed deconvolution technique for estimating IRAC 
fluxes in the ECDFS and examine the performance. 
The properties of the combined TENIS and IRAC 
catalog are discussed in Section [5] In Section [HI we 
summarize our results. The images and catalog from 
this work are availabl e at the official TENIS w ebsite: 
http : //www . asiaa. sinica. edu. tw/~bchsieh/TENIS/. 
All fluxes in this paper are /„. All magnitudes are in AB 
system, unless noted otherwise, where an AB magnitude 
is defined as AB = 23.9 - 2.51og(^Jy). 

2. CFHT WIRCAM IMAGING OBSERVATIONS 

The TENIS data were taken using WIRCam on the 
CFHT. WIRCam consists of four 2048 x 2048 HAWAII2- 
RG detectors covering a field-of-view (FOV) of 20' x 20' 
with a 0'.'3 pixel scale. The gap between the four detec- 
tors has a width of 45". The exposures were dithered 
to recover the detector gap as well as bad pixels. We 
adopted the standard dither pattern provided by the 
CFHT, which distributes the exposures along a ring with 
a default radius of 1/5. However, we changed the ra- 
dius every half semester or so, to 0.5 x to 1.5 x of its 
default value. This further randomizes the dither foot- 
prints and minimizes artifacts caused by flat fielding, sky 
subtraction, or crosstalk removal. At each dither point, 



we typically co-added two 50-second exposures and four 
20-second exposures at J and Ks, respectively. The typ- 
ical number of dither points were 9 for J and 11 for Ks, 
in a dither sequence. Therefore, after including the read- 
out time (10 seconds), a dither sequence took around 20 
minutes to complete. In our experience, the variation in 
sky color within such a short period is usually sufficiently 
small to allow for good flat fielding and sky subtraction. 

The ECDFS has a half-degree size, which is larger 
than the - 23'-25' dithered WIRCam FOV. To cover 
the ECDFS, we thus further offset the pointings between 
each dither sequence by ±2.'5 along RA and Dec. In the 
J-band imaging, we offset primarily along the NE-SW 
direction, leaving shallower NW and SE corners. In the 
ifs-band imaging, we were able to cover all the four cor- 
ners of the ECDFS. 

The observations were carried out by the CFHT in 
queue mode with weather monitoring. Most of the data 
were taken under photometric conditions with similar 
seeing. The average seeing for the observations are 0'.'8 
and Of! 7 (FWHM) for J and Ks, respectively. For the 
J-band data, we obtained 42.6 hours of on-source inte- 
gration in semesters 2007B and 2008B. For the -FCs-band 
data, we obtained 25.3 hours of integration in semesters 
2009B and 2010B. In addition, a Hawaiian group led by 
Lennox Cowie obtained 13.7 hours of Ks integration in 
semester 2010B. We include all the 39.0 hours of Ks data 
here. The integration maps for the J and Ks data are 
shown in Figure [21 and the cumulative areas in the TE- 
NIS J and Ks images as functions of effective integration 
time are shown in Figure [3) 

3. DATA REDUCTION 

The TENIS J and Ks data were processed using an In- 
teractive Data Language (IDL) based reduc tion pipeline. 
The d etails of the pipeline are described in lWang et al. I 
(2010), who also dealt with extremely deep WIRCam 
i maging. The reduc tion here is nearly identical to that 
in IWang et al. I pOlCft . 

We first grouped the images and only processed 
images from the same dither sequence and the same 
detector at a time. The dithered raw images were 
first median-combined to produce an initial flat field. 
Then objects were detected on each flattened image 
and were masked on the associated raw image. For an 
image, a better flat-field image was then created by 
median-combining the rests of the object-masked raw 
images in the dither sequence. This way, every image 
has an associated flat field, created by avoiding itself. 
This is the major difference between t he reduction 
here a nd that in the early data rel ease of IWang et al. I 
(|2010ft . In the earlier versions of the IWang et al. I (|2010T ) 
pipeline, all images in a dither sequence were median 
combined to form a flat field image, after detected 
objects were masked. This produces an extremely small 
but statistically significant overestimate of the flat field 
at the locations of faint objects that are undetected 
in single exposures. Because the sky background is 
very bright, this small error in flat field translates to a 
larger error in the final fluxes of faint objects. This is 
discussed in details in the later data release of Wang et al. 
(http : //www. asiaa. sinica. edu. tw/~whwang/goodsn_ks). 

6 Also see http: / /www. asiaa. sinica. edu. tw/~ whwang/idl/SIMPLE. 
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Fig. 1. — The TENIS J+Kg image is in grayscale. The green dot-dashed line indicates the ECDFS field, the red dotted line indicates 
the field of the SIMPLE data, and the blue dashed line indicates the field of the GOODS-S IRAC data. 
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Fig. 2. — Integration maps for the TENIS data. The left panel is for J and the right panel is for Kg. The grayscale shows the effective 
integration time in each panel where darker indicates a longer integration time. The white contours indicate 20%, 40%, 60%, and 80% of 
the maximum effective integration time, which is 42.6 hours for J and 39.0 hours for Kg. 



4 



Hsieh et al. 




Normalized Effective Integration Time 

Fig. 3. — Cumulative area in the TENIS J and K$ images vs. 
effective integration time. The maximum effective integration time 
(42.6 and 39.0 hours for J and Kg, respectively) is normalized to 
be 1.0 for each band. 

We would like to draw to the community's attention 
that a similar error may exist in other datasets if a 
similar flat-fielding method is adopted. Our reduction 
here avoids the above problem. 

After the above flat-fielding, usually the images are suf- 
ficiently flat and only a constant background subtraction 
is needed. Occasional residual sky structure caused by 
rapidly changing sky color was further subtracted by fit- 
ting fifth-degree polynomial surfaces to the image back- 
ground. 

On each WIRCam detector, there is crosstalk among 
the 32 readout channels (2048 x 64 pixels each). The 
crosstalk has different strength within the entire detector 
(32 channels), and within each of the four video boards 
(eight channels each). For every flattened and back- 
ground subtracted image, we removed the 32-channel 
crosstalk by subtracting the median combination of the 
32 2048 x 64 object-masked stripes. A similar procedure 
is then repeated to remove the 8-channel crosstalk. This 
greatly suppresses the crosstalk in the images, and only 
very weak residual effect can be seen in the final ultra- 
deep stack around several tens of the brightest objects in 
the ECDFS. 

To correctly stack the dithered images, optical distor- 
tion needs to be corrected. O ur reduction pipeline fol- 
lowed the method developed by I Anderson fc King] (2003) 
to derive the distortion function. The pipeline first de- 
tected objects in all dithered images, and then calculated 
the spatial displacement caused by the dithering for each 
object. Such displacement as a function of position in the 
images is actually the first-order derivative of the optical 
distortion function. The pipeline approximated the dis- 
placement function with polynomial functions of X and 
Y, and integrated them back to obtain the distortion 
function. We referred to Wang et al. (2010) for more 
detailed discussion about this technique. 

To obtain absolute astrometry and to project the im- 
ages onto the sky, we compared the distortion-corrected 
object positions with the source catalog produced by 
the HST Galaxy Evolution from Morphology and SEDs 



(GEMS, IRix et al. I l200i ICaldwell et al~l 12001 survey. 
By forcing the object positions to match those in the 
GEMS catalog, we computed a single two-dimensional, 
third-degree transformation function that contains all 
the effects including the distortion, sky projection, and 
absolution astrometry, for each dithered image. There- 
fore, in the entire reduction, each image only underwent 
one geometric transform. This minimizes the impact of 
image smearing caused by the transform. The trans- 
formed images can then be stacked to form a deep, as- 
trometrically correct image. 

Before images were stacked, photometry is carried out 
on individual exposures, and compared among the expo- 
sures. This allows to adjust relative zero-points of the 
individual exposures. After the images in a dither se- 
quence were stacked, absolute zero-point calibration was 
made by matching the fluxes in apertures of 5" in diam- 
eter with the "default magnitudes" in the point source 
catalog of the Two M icron All Sky Survey 6 (2MASS, 
ISkrutskie et all 120061 Q Since the magnitudes in the 
2MASS catalog are in the Vega system, we converted the 
2MASS magnitudes to fluxes based on the Vega zero- 
magnitudes in J and K$ of 1594 Jy and 666.8 Jy, re - 
spectively, as provided by 2MASS (jCohen et al. 112003ft . 
All the stacked, calibrated images from single-dither se- 
quences and from the four WIRCam detectors were then 
mosaicked for form a wide-field, ultradeep image. 

4. REDUCTION QUALITY 

4.1. Astrometry 

The GEMS survey used the Advanced Camera for Sur- 
veys (ACS) on the Hubble Space Telescope (HST) to 
image nearly the entire ECDFS. As described in the pre- 
vious section, we forced our astrometry to match that of 
the GEMS ACS catalog. Here we compare how well our 
astrometry matches that of the GEMS catalog. More 
than 9,000 compact sources with good detections (S/N 
> 20 and FWHM < 1'.'2) are selected for this comparison. 
Figures |4] and [5] show the relative astrometric offsets be- 
tween the TENIS J+K s image and the GEMS ACS cat- 
alog. As shown in Figure|H the systematic offsets in R.A. 
and Decl. are 0"016 and 0'.'0035, respectively. These off- 
sets are negligible compared with the rms scatter, which 
are 0'.'087 and 0'.'080 for R.A. and Decl., respectively. On 
the other hand, there are position-dependent systematic 
offsets. In the bottom- left panel of Figure El there is a 
wiggling structure with a scale of ~ 0.05 degree, which is 
compar able to the FOV of the ACS Wide-field Camera 
(WFC. lFordlt~al~ll200l FOV = 202" x 202"). Similar 
but less obvious structures can also be found in the other 
panels of Figure [5l It is highly likely that these position- 
dependent systematic offsets are internal to the GEMS 
catalog. 

According to Caldwel l et al. I (|2008l ). the astrometry of 
GEMS catalog is registered to that of the Classifying Ob- 
jects through Medium-Band Observat ions, a spectropho - 
tometric 17-filter survey rCOMBO-17. lWolf et al. 1l200l . 
The absolute astrometric accuracy of the GEMS cata- 
log is therefore limited by the astrometric quality of the 

7 2MASS is a joint project of the University of Massachusetts 
and the Infrared Processing and Analysis Center/California Insti- 
tute of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 
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Fig. 4. — Relative astrometric offsets between the TENIS J+Ks 
image and the GEMS ACS catalog. We selected over 9,000 compact 
sources with good detections (i.e., S/N > 20 and FWHM < 1'.'2) 
for the comparison. The systematic offsets in R.A. and Decl. are 
CK'0016 and 0"0035, which are negligible as compared to the rms 
offsets, which are 0'.'087 and / /080 for R.A. and Decl., respectively. 

COMBO-17 catalog, which is better than 0'.'15 but which 
may be greater than 0'.'3 in some localized regions. As 
we calibrated our astrometry using the GEMS catalog, 
the absolute astrometric accuracy of the TENIS data is 
therefore also limited by that of the COMBO-17 data. 
We could match the TENIS astrometry directly to the 
COMBO-17 data. However, the major scientific goal of 
the TENIS project is to find LBGs at z > 7, which needs 
the ultra deep GEMS data. We therefore matched the 
TENIS astrometry to that of the GEMS data rather than 
that of the COMBO-17 data because it makes combining 
the TENIS data and the GEMS data much easier. 

4.2. Photometry 

As mentioned in Section El the TENIS WIRCam J and 
Ks fluxes were calibrated using 5" diameter apertures 
and matched to the fluxes in the 2MASS point source 
catalog. In Figure HI we show the flux ratios between 
2MASS and TENIS in J and K s , where the flux range 
between the two vertical dashed lines in each panel indi- 
cates that used for the calibration. Fluxes brighter than 
this range suffer from nonlinearity issues in WIRCam, 
whereas fluxes fainter than this range are subjected to 
selection effects due to the much shallower detection lim- 
its of 2MASS. The error- weighted means of the objects in 
the chosen calibration flux ranges are 1.0002±0.003 and 
1.0006±0.005 for J and Ks, respectively, resulting in an 
overall flux calibration quality that is good to 0.3% and 
0.5% for J and Ks, respectively. We note that the flux 
errors of the TENIS data are also shown horizontally in 
Figure [6] but are too small to be visible, and hence they 
make negligible contributions to the photometric error 
budget compared with the vertical error bars. The flux 
calibration quality of the TENIS data is therefore com- 
pletely dominated by the 2MASS flux errors. 

We also investigated the photometric uniformity in our 
TENIS images by subdividing the TENIS J and K s im- 
ages into four quadrants and then comparing the TE- 



TABLE 1 

SExtractor Parameters 



Parameter 


Value 


DE 1 EC 1 JV11N AREA 


2 


DEI EC 1 _ 1 HRESH 


1.3 


A TVT A T VCTQ HPUTD T7GU 


1.3 


b 1L1 hjH 


Y 


T7TT 'T TT T'> TVT A ft JT? 

b 1L 1 hjH-lN ANLbj 


gauss_2 . 0_3x3 . conv 


DEBLEND.NTHRESH 


64 


DEBLEND.MINCONT 


0.00001 


CLEAN 


Y 


CLEAN.PARAM 


0.1 


BACK.SIZE 


24 


BACK.FILTERSIZE 


3 


BACK.TYPE 


AUTO 


BACKPHOTO.TYPE 


LOCAL 


BACKPHOTO.THICK 


40 


WEIGHT.TYPE 


MAP.WEIGHT 


WEIGHT.THRESH 


20 



NIS fluxes with the 2MASS fluxes as we did for the en- 
tire images. For quadrants 1 through 4 in J, the error- 
weighted flux ratios are 0.979 ± 0.007 (21 sources), 0.997 
± 0.006 (25 sources), 1.014 ± 0.005 (43 sources), and 
1.005 ± 0.008 (22 sources), respectively. For those in 
K s , they are 0.998 ± 0.009 (19 sources), 1.033 ± 0.007 
(27 sources), 1.003 ± 0.009 (27 sources), 0.978 ± 0.011 
(11 sources), respectively. The standard deviations of the 
four measured offsets are 1.3% and 2.0% for J and Ks, 
respectively. This suggests that the photometric gradi- 
ents over size scales of ~ 15' are less than 0.013 and 0.02 
mag in the TENIS J and Ks images, respectively; these 
measured 1.3% and 2.0% are fairly good as compared 
to other NIR extragalactic deep surveys ov er similar size 
scale s (e.g., 2% in the COSMOS survey; iCapak et al. I 
[2001 . 

5. THE J AND K s PHOTOMETRIC CATALOG 

To produce a source catalog that is as complete as pos- 
sible, we decide to perform object detection in a high S/N 
image generated by combining the TENIS J and Ks im- 
ages. Such a combination has to take into account the 
different integration time distributions in J and Ks- A 
straightforward method is to weight the pixels by their 
associated integration times and then combine the im- 
ages. However, our J integration time is approximately 
1.32x longer than the Ks integration time, but objects 
are generally brighter (in our map unit, which is /^Jy) 
at Ks- Hence, a direct integration time weighted com- 
bination is less optimal in terms of combined S/N for 
most objects. We therefore normalized the J and Ks in- 
tegration times by artificially reducing the J integration 
time by a factor of 1.32 (i.e., downweighting the J im- 
age), and then performed the integration time weighted 
J+Ks combination. 

We used SExtractor version 2.5.0 (Be rtin fc Arnoutsl 
1996) to detect objects and measure their fluxes in the 
TENIS WIRCam images. The "FLUX.AUTO" values 
of the SExtractor output are chosen for the flux mea- 
surements. The double-image mode of SExtractor was 
performed to detect objects in the J+Ks image while 
measuring fluxes in the original J and Ks images at 
the locations of the J+Ks objects. The most important 
SExtractor parameters we used are listed in Table [1] 

The flux errors provided by SExtractor are derived 
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procedure. First, the fluxes and the flux errors were re- 
measured using SExtractor with 2" diameter apertures. 
We then convolved the source-masked image with a 2" 
diameter circular top-hat kernel, and calculated the rms 
in a 10" x 10" area around each pixel on the convolved 
image. The ratio between the aperture photometric flux 
error provided by SExtractor for a certain object and 
the rms value around the same position in the convolved 
image is the correction factor for its flux error. The me- 
dian value of the correction factors was computed to be 
the general correction factor for all the sources. For the 
TENIS J- and ifg-band data, the general correction fac- 
tors are very similar, which is 1.27. The same procedure 
was repeated with different aperture sizes. We found 
that the correction factors are very stable over different 
ap erture sizes, which is consistent with the experience 
in iWang et al. I (|2010[) . Therefore we applied this fac- 
tor to the "FLUXERR_AUTO" values for all the objects 
to account for the effects of noise correlation and confu- 
sion. We note that this factor does not apply to the flux 
measurements. 

We then estimated the aperture correction for the 
AUTO aperture. Since wc calibrated the WIRCam data 
to 2MASS using a 5" diameter aperture, we just need to 
derive the correction factor between the AUTO aperture 
and the 5" aperture. After comparing the fluxes derived 
using these two different apertures, we found that they 
are in very good agreement with each other for most 
objects. Because the photometry derived using AUTO 
apertures have better S/N for objects with various differ- 
ent morphologies and fluxes as compared to the 5" fluxes, 
we decided to use the original "FLUX_AUTO" values as 
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Fig. 6.— Flux ratios between 2MASS and TENIS in J and K S - 
The two vertical dashed lines in each panel indicate the flux range 
for the flux calibration. We note that the outliers between the two 
vertical dashed lines are excluded in the flux calibration. The flux 
errors of the TENIS data are also shown horizontally but they are 
too small to be visible. Hence they have negligible contribution to 
the vertical error bars. The flux calibration quality of the TENIS 
data is therefore completely dominated by the 2MASS flux errors. 



from the background noise directly, and the noise co- 
variance between pixels (i.e., correlated error produced 
by image resampling, and to a lesser degree, faint unde- 
tected objects) is ignored. To mitigate against the lat- 
ter effect, we calibrated the errors using the following 
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the final flux measurements in our catalog. 

6. COMPARISON WITH THE GOODS-S/ISAAC DATA 

There are deep near-infrared imaging observations of 
the GOODS-S region using the Infrared Spectrometer 
And Array Camera (ISAAC) on the Very Large Tele- 
scope (VLT) in J, H, and Ks- The jig -selected cata- 
log is published by iRetzlaff et al. I (120101) . Accordi ng to 
the sensitivities provided by Rctzla ff et al. I (|2010h . the 
depths of the ISAAC data are comparable with that of 
TENIS. We have checked the photometric consistency 
and the data quality difference between the two catalogs. 
In Figure [71 we compare the TENIS J and Ks photome- 
try with the total magnitudes in the ISAAC catalog. The 
ISAAC fluxes for bright stars are approximately 10% and 
15% less than the TENIS fluxes at J and Ks, respec- 
tively. We checked whether these large differences could 
be caused by the differences in the filter systems used, 
as shown in Figure [5J We found that the different pass- 
bands of the two filter systems can produce only ~ 3% 
differences in fluxes by ex amining the co lor of stars from 
Kurucz models (ATLAS9; lKurucz 1119931 ). It is also worth 
noting that both catalogs are not corrected for Galac- 
tic extinction, although the correction values, only 0.008 
mag (J) and 0.003 mag (Ks), are too small to explain 
the > 10% flux differences. 

We calibrated the TENIS photometry using the 
2MASS catalog. The ISAAC observations, however, were 
calibrated with standard stars. We therefore directly 
compared the ISAAC catalog with the 2MASS catalog 
to see if the flux differences are due to different zero- 
point calibration methods used in the two catalogs. The 
result is shown in Figure |H The ISAAC total fluxes are 
~ 10% to 15% lower than the 2MASS default fluxes, 
consistent with the differen ces between the TENI S and 
ISAAC fluxes. In addition, IRetzlaff et al~l ([20101 ) men- 
tioned that a significant bias (<0.1mag) is visible when 
comparing ISAAC to the 2MASS catalog. Hence we con- 
clude that the flux discrepancies between the TENIS and 
ISAAC catalogs are due to different zeropoint calibration 
methods. 

To compare the depth, we plot the magnitude vs. S/N 
for the TENIS and ISAAC data in Figure [TO] The 5a 
limiting m agnitudes for point so urces of the ISAAC data 
claimed by IRetzlaff et al. I ()2010D are about 25.0 and 24.4 
for J and Ks, respectively, which are consistent with 
what are shown in Figure 1101 According to this fig- 
ure, the TENIS data are about 0.5 mag deeper than 
the ISAAC data in both J and Ks- If the abovemen- 
tioned zeropoint differences are considered, the depth 
differences would be further increased to > 0.6 mag. 
We note that the Cosmic Assembly Near- infrared Deep 
Extragalactic Legacy Survey (CANDELS, iGrogin et al. I 
120111 : iKoekemoer et al.ll2011l ) provides a deeper J-band 
(F125W) dataset in the GOODS-S region (5a limiting 
magnitudes for point sources is ~ 27 mag). We therefore 
conclude that our TENIS catalog is by far not only the 
deepest NIR dataset in ECDFS, but also the deepest Ks 
dataset even in the narrow GOODS-S region. 

7. IRAC PHOTOMETRY 

The SIMPLE survey (|Damen et al. I 120111 ) provides 
deep IRAC observations covering the entire ECDFS with 
the 10' x 15' ultradeep GOODS-S IRAC (M. Dickinson et 



al., 20 12, in preparation) m osaics in its center. Accord- 
ing to iDamen et al. I ([201 ll ) , the 5a limiting magnitudes 
for point sources typically are 23.8, 23.6, 21.9, and 21.7 
for 3.6/zm, 4.5/xm, 5.8/^m, and 8.0/^m, respectively. For 
the ultradeep GOODS-S IRAC region, the 5a limiting 
magnitudes for point sources are 26.1, 25.5, 23.5, and 
23.4 for 3.6/ zm, 4.5/zm, 5.8 A tm, a nd 8.0/im, respectively, 
according to Dahl en et al. I (|2010[) . As mentioned in Sec- 
tion [ll IRAC data are valuable for studies from Galac- 
tic objects to the high-redshift universe. Combining the 
SIMPLE IRAC and TENIS data can broaden the use of 
the TENIS catalog. 

7.1. Basic Principle 

The major difficulty in combining the IRAC and 
WIRCam data is the relatively large IRAC PSFs (T/5 
to 2'.'0) compared with of the WIRCam data (0'.'8). As 
a consequence, measurements of total fluxes of objects 
in the IRAC images require photometric apertures that 
are much larger than the WIRCam apertures. Given the 
ultradeep nature of the SIMPLE and GOODS-S IRAC 
images and the high surface densities of objects, large 
apertures mean that the photometric accuracy is highly 
subject to nearby bright objects as well as the large num- 
ber of faint confusing sources. Because the blended ob- 
jects can have very different morphology and color, ac- 
counting for this problem is crucial. 

Many intensive efforts have gone into better esti- 
mating fluxes in crowded, low-resolution images. The 
most recent methods rely on utilizing the positional 
and morphological information of objects detected in 
a high resolution image in a different waveband (e.g. 



lLaidler et al. I l2007t iWang et all 



Gra zian et al. 1 120061 

120101: iMcLure et al. l[20rl . By assuming that the in- 
trinsic morphology of an object is identical in the two 
wavebands (i.e., no color gradient from center to edge) 
and taking into account differences in the PSFs, one can 
model how an object would look like in the low reso- 
lution image based on its observed morphology in the 
high resolution image . Usi ng a minimum y 2 m ethod, 
iGrazian et"al~l ([20061 ) and lLaidler et al~l |2007) were 
able to model all blended objects simultaneously and 
thus estimate their fluxes in the low resolution image. 
Instead of minimizing y 2 , another approach is to mini- 
mize residual fluxes be tween the model an d the observed 
low resolution image ([Wang et al. 1 120101 ) . This is the 
basic idea of CLEAN deconvolution in radio imaging 
([Hogbom et al. Ill97l . 

In this w o rk, w e follow the CLEAN approach of 
IWang et al. I ([20101 ). However, we greatly relax the as - 
sumption on object morphology. In Wa ng et al. I (|2010f >. 
a model of the IRAC image of a galaxy is generated based 
on the PSFs and its high-resolution WIRCam image, and 
the model is then iteratively subtracted ("CLEANed") 
from the real IRAC image. Here, we directly CLEAN 
the IRAC image of a WIRCam detected galaxy using 
the IRAC PSF, and we allow the cleaning position to 
move around its WIRCam position. In other words, our 
method is essentially the same as CLEAN deconvolution 
in radio imaging, with nearly no restrictions except for 
the locations where CLEAN can happen. We named this 
method "IRACLEAN" because it is designed for estimat- 
ing IRAC fluxes. 
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Fig. 7. — Photometry comparisons between the TENIS and ISAAC catalogs. Flux ratios between the two catalogs versus the TENIS 
fluxes are plotted. The left panel is for J and the right panel is for Kg. The red line indicates the running median and the green lines 
indicate the upper and lower 68th percentiles. The blue dashed line indicates a -0.15 mag difference. The results show that the ISAAC 
fluxes for bright stars are about 10% to 15% less than the TENIS fluxes. 
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Fig. 8. — Transmission curves for the WIRCam and ISAAC filters 
is for WIRCam and the dashed curve is for ISAAC in each panel. 

7.2. Methodology 
7.2.1. PSF Construction 

The IRAC PSFs were generated using bright isolated 
point sources in each IRAC channel. The size of the 
PSF images is V x V in order to cover the outer struc- 
ture of the extended IRAC PSFs, especially for 8.0/im. 
To avoid saturated sources, we chose intermediate bright- 
ness objects and checked their images carefully. We also 
checked their optical counterparts in the extremely high 
resolution HST ACS images to make sure that they are 
really point-like. The final PSF is constructed by cal- 
culating the flux weighted mean of the PSFs of these 
objects with a 3ct clip. There are, however, still many 
objects within 1' x 1' around these sources even if they 
are relatively "isolated" . This biases the wing of the con- 
structed PSF. We therefore ran SExtractor on the IRAC 
images to generate object masks and then masked all the 
nearby objects when generating the PSFs. This masking 



E 




2.1 2.2 2.3 2.4 
Wavelength (urn) 

The left panel is for J and the right panel is for Kg. The solid curve 

procedure, however, may create many "holes" in the con- 
structed PSF image. To overcome this, sometimes more 
than 20 objects have to be used for constructing one PSF. 
The PSFs of stars in the native IRAC images would be 
under-sampled because of the relatively large pixel size of 
IRAC (l'/2 per pixel for the IRAC PSFs with FWHMs of 
l'/66 to l'/98). Without doing a sub- pixel centering, the 
constructed PSF may be artificially broadened. Since the 
SIMPLE image has been resampled to match the TENIS 
image (~ 0'.'3 per pixel), simply finding the location of a 
peak is very similar to do a sub-pixel centering, and the 
centering accuracy of the constructed PSF is better than 
9% to 7.6% (±0'.'15, from 3.6/xm to 8.0/mi) of the IRAC 
FWHMs. Hence the broadening effect of the stacked PSF 
is negligible. We note that the PSFs in the SIMPLE im- 
ages are position-dependent, and the procedure t hat we 
used to deal with this issue is described in Section [7.2.5! 

7.2.2. Object Identification 
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Fig. 9. — Same as Figure but for 2MASS and ISAAC. The blue dashed line indicates a -0.15 mag difference. The flux differences 
are consistent with that shown in Figure [7] which suggests that the flux discrepancies between the TENIS and ISAAC catalog are due to 
different zeropoint calibration methods. 
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Fig. 10. — S/N vs. Magnitude plots for the TENIS and ISAAC data. The left panel is for J and the right panel is for Ks- The dark 
black dots indicate objects in the GOODS-S region in the TENIS catalog and the gray dots indicate objects in the ISAAC catalog. The 
horizontal line in each panel marks the S/N level of 5. Based on this comparison, the TENIS data are about 0.5mag deeper than the 
ISAAC data in both J and Ks- 



We used the TENIS J+Ks image as a prior to esti- 
mate the IRAC fluxes around the locations of J+K$ de- 
tected objects. The location information is contained in 
the SExtractor "segmentation map" of the TENIS J+K$ 
image. The segmentation map tells which detected ob- 
ject (or no object) a pixel is associated with. For a given 
J+Ks detected object, we only clean its IRAC fluxes 
at its associated segmentation map pixels. Therefore, 
the boundary of an object defined in the segmentation 
map is equivalent to "CLEAN window" in radio imag- 
ing. In addition, since we allow any pixels within the 
boundary of an object to be CLEANed, we do not as- 
sume any morphology for that object except for its outer 
extent as defined by the J+Ks image. This has the ad- 



vantage of allowing changes in morphology at different 
wavelengths. We should note that changing the SEx- 
tractor detection threshold setting for the J+Kg image 
(i.e., DETECT.THRESH) would enlarge and shrink the 
segmentation area, which might affect the IRACLEAN 
performance. To ensure IRACLEAN working properly, 
the detection threshold should be set as low as possible. 
Lowering the detection threshold down too much, how- 
ever, would also increase the spurious rate significantly. 
We found that the IRACLEAN fluxes with different DE- 
TECT.THRESH values agree with each other within 3% 
as long as DETECT.THRESH < 1.5, and the spurious 
rate increases dramatically with DETECT.THRESH < 
1.2 (see Section [8T21 for details about spurious rate). Set- 
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ting DETECT.THRESH = 1.3 (see Tabled]) is a good 
balance between IRACLEAN performance, source com- 
pleteness, and spurious rate. 

7.2.3. IRACLEANing Objects 

We first made mosaics of the IRAC images and resam- 
pled them to match the pixel size of the TENIS J+Ks 
image (0'.'3). The IRACLEAN process always starts at 
an IRAC pixel with the highest "absolute" value (i.e., the 
value can be negative) measured within a 9 x 9-pixel box 
(Fap9, hereafter), and this pixel (Pdecon, hereafter) must 
have been registered to one object in the J+Ks segmen- 
tation map. We chose a 9 x 9-pixel (2'.'7 x 2'.'7) box as the 
aperture size because it delivers the best S/N. Since the 
IRAC image (0'.'6 per pixel) has been resampled to match 
the TENIS J+Ks image (0'.'3 per pixel), simply moving 
a window (i.e., F^pg) across the IRAC image to find the 
location of a peak is very similar to doing a sub-pixel 
centering. Once Pdecon is found, we subtracted a scaled 
PSF from the surrounding 1' x 1' area centered on Pdecon- 
The percentage of the subtracted flux ("CLEAN gain") 
depends on Fapq, if Fapq is greater than 20er@, then 1% 
of the flux was subtracted; if Fap9 is less than 20cr but 
greater than 5a, then 10% of the flux was subtracted; 
if Fap9 is less than 5a, then 100% of the flux was sub- 
tracted. Our gain for brighter peaks (1% for 20cr) is much 
smaller than the values in most radio CLEAN (10%). 
We found a gain of 10% for saturated or extended bright 
sources sometimes produces unreliable results. This is 
likely because the IRAC PSF is not as well determined 
as the synthesized beam in radio imaging, so that over- 
or under-subtraction easily occurs. We found that a gain 
of 1% produces a good balance between processing speed 
and deconvolution quality. The subtracted flux (Fsub) 
was summed and registered to the associated object in 
the J+Ks segmentation map. After each subtraction, 
the IRACLEAN process is repeated on the subtracted 
image, until there were no pixels with Fap9 higher than 
2a. This threshold is chosen such that it is not too low 
to allow for unreliable objects entering IRACLEAN, and 
it is not too high to cut off too much useful information 
on weak objects. We emphasize again that the concept 
of the deconvolution here is identical to that of CLEAN 
in radio imaging. 

Under some extreme conditions (unmatched PSF with 
saturated or bright extended objects), gains of < 1% 
still do not work. The residual fluxes of these objects 
would start to oscillate and diverge with the progress 
of IRACLEAN, and this phenomenon usually happens 
when the residual fluxes are about 30% to 50% of the 
original fluxes. Since the PSF area we used is 1' x 1', this 
effect would also seriously affect the IRACLEAN results 
of other objects within the 1' area. The IRACLEAN 
process will set a gain of 100% to subtract the fluxes 
for such objects when the signs of flux divergence show 
up, and then stop to clean them hereafter. It is worth 

8 The sigma here is the local background fluctuation measured 
using the following steps: (a) running SExtractor in the single- 
image mode on all the IRAC mosaics and generating the segmen- 
tation maps; (b) masking the detected sources according to the 
segmentation maps to generate the background images; (c) con- 
volving the background images with a 9 X 9-pixel top-hat kernel 
and then generating the noise maps by calculating the rms around 
each pixel on the convolved background images. 



noting that since greater than 50% of the flux is still 
cleaned by the normal IRACLEAN procedure, the bias 
effect of assigning a gain of 100% to the residual fluxes 
should be much less than 50% (i.e., <C0.5 mag). 

7.2.4. Flux and Error Measurements 

At the end of the IRACLEAN process, the summation 
of Fsub for each object is the flux measurement of that 
object, and the final subtracted image is the residual 
map. If an object cannot be detected in IRAC images 
(i.e., its Fap9 is lower than 2a), 0.0 is assigned to its flux 
measurement. The residual map allows us to check the 
quality of IRACLEAN, and is also used for estimating 
flux errors. The flux error of each object was calculated 
based on the fluctuations in the local area around that 
object in the residual map. Any imperfection of the PSF 
would cause larger fluctuations in the residual map, and 
this effect is included in the flux error calculation. 

To estimate the flux error of an object, we convolved 
the residual map with a 9 x 9-pixel top-hat kernel and 
then generated a noise map by calculating the rms 
around each pixel on the convolved residual map. In con- 
ventional aperture photometry, the flux error scales with 
the square-root of the aperture area. To follow this rule, 
we thus defined an effective aperture size for each object 
using the following method. First, if the IRACLEAN 
process always deconvolves an object at the same Pdecon 
(i.e., deconvolving a bright point source with a matched 
PSF), then the effective aperture size for this object is 
81 pixels (9 x 9, FigureEJi). If the IRACLEAN process 
deconvolves a slightly extended object on two adjeent 
Pdecon-, then the effective aperture size for this object can 
be 90 pixels (Figure ITTb ) , or 98 pixels (Figure \TTb ) . For 
a more extended source, the distribution of these Pdecon 
can be more complicated, and the shape of the aperture 
can be irregular ( Figure [TTtO ■ 

Once the effective aperture size is determined, the flux 
error can be estimated with 

n 

AaT — ^ rmSl 

where N p i xe i is the total number of pixels in the effective 
aperture, n is the total number of Pdecon, and rmsi is the 
value from the noise map for a certain Pdecon- The first 
term in Equation [T] is to scale the flux error with the ef- 
fective aperture size, and the second term is the average 
flux error for the local area around this object. One can 
use a more complicated method of calculating the sec- 
ond term by compute the local flux errors using different 
weightings. For simplicity, we just use an unweighted 
mean, which is similar to that used in conventional aper- 
ture photometry (cf., e.g., photometry with PSF fitting). 
Using Equation [TJ we are able to reasonably estimate the 
flux error for all objects regardless of their morphologies. 
For an object undetected in an IRAC image, we used the 
coordinate of that object derived from the J+Ks image, 
and calculated its flux error as if it was a point source 
(i.e., Using Equation [T] with N p i xe i = 81, n = 1, rms 
= the value from the noise map for the Pdecon, where 
the Pdecon is the object coordinate in the J+Ks image 
provided by SExtractor). 

Since we used a 9 x 9-pixel box as the flux aperture for 
IRACLEAN, an aperture correction needs to be applied 
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Fig. 11. — Examples for calculating aperture sizes for various 
cases. The grid in grayscale indicates the pixels in the resamplcd 
IRAC images; the pixel marked with black color indicates Pd econ ; 
and the area marked with red color shows the aperture size for this 
object. Four different cases are shown in this plot: (a) a bright 
point source (only one Pdecon i n the object center): the aperture 
shape is a 9 X 9-pixcl box, i.e., the total number of pixels in the 
aperture is 81. (b) a slightly extended source (two nearby P,j econ s 
around the object center): the aperture shape is a 9 X 10-pixel 
rectangular, and the total number of pixels in the aperture is 90. (c) 
a slightly extended source with different position angle as compared 
to case (b): the aperture size is 98 pixels, (d) an extended source: 
the aperture shape is more irregular and the aperture size is 140 
pixels. 

for the flux errors. The value of aperture correction is de- 
rived from the deconvolution PSF for each channel. We 
describe the aperture correction values in Section 17.2.51 

7.2.5. Position-Dependent PSF 

The PSFs in the ECDFS IRAC images change with 
position. This is mainly caused by the change in the 
orientation of Spitzer during the observations. In the 
GOODS-S area, the upper and lower half fields were ob- 
served with IRAC in two different epochs, with some 
overlap in the middle. The orientation of Spitzer 
changed by nearly 180 degrees between the two epochs. 
The SIMPLE data were also taken in two different epochs 
with a several-degree difference in Spitzer's orientation. 
More than 80% of the SIMPLE area is covered by both 
epochs. In addition, there is also some overlap between 
the GOODS-S and SIMPLE observations. Because of 
the asymmetric IRAC PSFs and different orientations 
of Spitzer in the various observing epochs, the stacked 
IRAC image has several different synthesis PSFs. Using a 
single PSF to clean the entire image will significantly de- 
grade the IRACLEAN quality because the performance 
of IRACLEAN critically depends on the accuracy of the 
deconvolution PSF. IRACLEAN is affected by the accu- 
racy of the deconvolution PSF in two aspects: (a) For a 
bright point source, using an inaccurate PSF would lead 
to distorted intermediate residual images after iterating 
the clean process many times. From this point, the lo- 
cal brightest pixel (i.e., Pdecon) may not be around the 
center of this object so that IRACLEAN starts treat- 
ing it as an extended source until the end of the clean 
process. This effect would lead to a biased flux mea- 
surement. Same thing would happen to bright extended 
sources, too. Fainter sources are not affected by this ef- 



fect because their corresponding segmentation areas are 
only a few pixels hence their Pdecon pixels cannot go too 
far from the object centers, (b) Some fainter sources are, 
however, affected by another issue because of using an 
inaccurate PSF. If they are close to another objects that 
are not cleaned using accurate PSFs, then the dirty resid- 
uals of the wings of their close neighbors would contan- 
minate their flux measurements. The brighter the close 
neighbor is, the more serious the effect is; the closer the 
separation is, the more serious the effect is. In order to 
obtain better flux measurements for all objects, we used 
several different PSFs to clean one IRAC image. 

It is possible to construct many different PSFs in an 
IRAC image. We start with as many as seven different 
PSFs but below we will demonstrate that only four PSFs 
are necessary for 3.6 and 4.5 /an. First, we generated 
seven PSFs by picking up bright point-like sources from 
seven different areas in the IRAC image for each IRAC 
band. The generated PSFs are shown in Figurc[T2J Then 
for a given IRAC band, we repeated the IRACLEAN pro- 
cess seven times using one of the seven PSFs each time. 
This gives seven flux measurements for each object. Since 
our flux errors (Section 17.2.41) contain the effect of un- 
matched PSF, we can use the flux errors as an indicator 
on which PSF works the best on each object. 

Figures [TS] and [TJ] compare the flux errors based on 
the seven PSFs at 3.6 and 4.5 /xm, respectively. In the 
z-th column, we pick up objects whose flux errors are the 
lowest when IRACLEANed with the i-th PSF. Then in 
the j-th row, we compare their flux errors based on the 
j-th PSF with respect to the flux errors based on the i- 
th (best) PSF. The x-axis in each panel is object ID and 
the y-axis is the error ratio of j-th to i-th PSF. (Since 
in each panel we are always comparing against the best 
PSF, the flux ratios in all the panels are always greater 
than 1.) In other words, the comparison shows how the 
other PSFs perform when compared to the best-matched 
PSF for each object. 

By examining Figures [13] and 1141 we can determine 
which PSF is necessary and which PSF is redundant. In 
a panel, if most of the error ratios are fairly close to 1, it 
means that the IRACLEAN results using the i-th PSF 
and the j-th PSF have similar quality. In such a case, the 
i-th and j-th PSFs can be replaced by each other. For 
example, we found that for 3.6 /j,m, PSF5 and PSF7 can 
both be replaced by PSF2, and PSF1 can be replaced by 
any of the other PSFs. In other words, only PSF2, 3, 4, 
and 6 are necessary for IRACLEAN at 3.6 /im. Similarly, 
only PSF1, 2, 4, and 5 are necessary for 4.5 /xm. For 5.8 
and 8.0 /im, we do not find significant differences among 
the seven PSFs, so only one PSF each is adopted and is 
shown in Figure [T2l 

To show how seriously IRACLEAN is affected by the 
accuracy of the deconvolution PSF, we calculated the 
rms scatters of the IRACLEAN measurements using the 
seven different PSFs for all the objects in 3.6 /itm and 
4.5 /xm, and compared them with the noises measured 
from the residual images. The results are shown in Fig- 
ure[T51 According to FigureUHl the IRACLEAN fluxes of 
the bright isolated objects are affected by PSF more se- 
riously than that of the faint isolated objects, which can 
be explained by the abovementioned effect (a). Most of 
non-isolated objects are faint objects since bright objects 
have lower surface number density hence they are more 
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Fig. 12. — The PSF images for IRAC 3.6^tm to 8.0/xm. The image size is 1' X 1' in each panel. There are seven testing PSFs for 3.6/^m 
and for 4.5/Ltm, and one PSF for 5.8/im and 8.0/im. Contours are plotted on the PSF images for IRAC 3.6/^m and 4.5/rni to make the 
differences between PSFs visually clearer. 



isolated. These non-isolated faint objects show another 
distribution in this plot. Many of them are also seriously 
affected by inaccurate PSF, which can be explained by 
the abovementioned effect (b). Figure 14 proves that 
the performance of IRACLEAN critically depends on the 
PSF, and using multiple PSFs to clean the IRAC 3.6 /xm 
and 4.5 /im images can gain the most advantage from 
IRACLEAN. 

With the above experiments, we adopted the IRA- 
CLEAN results based on the final four PSFs for 3.6 and 
4.5 /xm, and one PSF each for 5.8 and 8.0 /xm. For a 
given object, we pick the flux measurement with the low- 
est error among the four PSFs at 3.6 and 4.5 /im. We 



then calculated the aperture correction factors for each 
channel based on the PSFs. For 3.6 /xm and 4.5 /xm, the 
correction factors derived from the four PSFs agree with 
each other within 1% and we adopted the mean values. It 
is worth noting that the correction factors from different 
PSFs are very similar because these PSFs are synthesized 
from the same PSF with different orientations, hence the 
ratios of fluxes between inside and outside the aperture 
for different PSFs agree with each other very well. It also 
suggests that the centering issue of generating the stack- 
ing PSFs is negligible. The adopted aperture correction 
factors are 1.894, 1.875, 2.111, and 2.334 for 3.6, 4.5, 5.8, 
and 8.0 /xm, respectively, for our 2'.'7 x 2'.'7 flux aperture. 
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Fig. 13. — Comparison of the performance of the seven PSFs at 3.6 ^m. The i-th column contains objects whose flux errors are the 
lowest when IRACLEANed with the i-th PSF. The j'-th row compares the errors derived with the j-th PSF against the errors with the i-th 
(best) PSF. X-axis is object ID and Y-axis is error ratio in each panel. See text for details. 



7.2.6. Known Issues 

Three extremely bright (thus saturated) stars and 
several very bright objects are in the ECDFS. They 
have very bright wings and occupy areas that are much 
broader than 1' x 1' (the deconvolution PSF size). Ob- 
jects close to these bright sources, including those around 
the diffraction spikes and crosstalk features, therefore 
have higher flux measurements. In addition, as men- 
tioned in Section 17.2. 31 we stop the IRACLEAN process 
if it starts diverging for some bright objects. This also 
implies that the flux measurements of nearby objects are 
affected. However, we should note that the flux measure- 
ments of these nearby objects will be affected even more 
seriously if we do not stop the IRACLEAN process when 
divergence happens. 

7.3. Quality and Performance 
7.3.1. Residual Images 



We demonstrate the performance of IRACLEAN in 
Figure El Regions of sizes 300" x 140" in the SIM- 
PLE IRAC images (left panels) and their residual images 
(right panels) are shown. All images are shown with in- 
verted linear scales. The brightness and contrast of each 
panel are identical. By visually checking the residual im- 
ages, we found that IRACLEAN works reasonably well in 
all four channels. For 3.6/zm and 4.5/im, however, there 
are halos around many sources in the residual images, 
which may be due to unmatched deconvolution PSFs 
and/or under-sampled PSF caused by the large pixel 
scale of IRAC. The 3.6 /nm or 4.5 /zm residual images 
shown in Figure [16] are generated using a single deconvo- 
lution PSF. However, as discussed in Section l7.2.5l there 
are indeed four PSFs for each band, and thus three more 
versions of residual images each, which are not shown 
here. Since we adopt the result that has the lowest resid- 
ual fluctuation for each object, the IRAC fluxes for 3.6/im 
and 4.5/mi in our photometric catalog are indeed better 



14 



4 
3 
2 
1 
4 
3 
2 
1 
4 
3 
2 

2 1 
1 4 

5a 

?2 

m 1 

4 
3 
2 
1 
4 
3 
2 
1 

4 
3 
2 
1 



Hsieh et al. 

PSF1 PSF2 PSF3 PSF4 PSF5 PSF6 PSF7 



i i i i i I i i i i i i i I i i i r 



1 1 1 1 1 1 1 



Min i 



1 1 1 ' ' 1 i I ' ' ■■ 1 1 ■ 



i l i l i l i 



M ill i i i i i 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



1 1 1 1 1 1 1 



"i i i l i i i r 



fffm 



Mm ! 



I l l l l l l 







D 

CO 



"0 

co 



n 

co 



TO 

co 



u 

co 



Object ID (4.5 (am) 



Fig. 14. — Same as Figure [13] but for 4.5/im. 

than the visual impression based on any single residual 
image. On the other hand, for 5.8/im and 8.0/mi, the 
residual maps are very clean; there are no obvious resid- 
ual effects or artifacts. 

7.3.2. Monte-Carlo simulations 

To further understand the performance of IRACLEAN, 
we carried out simple Monte-Carlo simulations for ob- 
jects with high S/N (S/N~ 50) blended by nearby ob- 
jects. All the simulated objects have flat spectra and 
are point-like. The separations between the blended ob- 
jects are from V.'2 to 3'.'9 and the flux ratios between 
the blended objects are from 1 to 100. We simulated a 
J+Ks image of the blended objects and corresponding 
images for four IRAC channels. The PSF of the J+Ks 
image is extracted from the real J+Ks data. The PSF 
for each object in the simulated 3.6/im and 4.5/im images 
is randomly picked up from the seven PSFs described in 
Section 17.2.51 and in Figure rT2J We ran SExtractor on 
the simulated J+Ks images to generate the segmenta- 
tion map, and then performed IRACLEAN on the sim- 



ulated IRAC images. We adopted only four PSFs for 
IRACLEAN for 3.6/im and for 4.5/im, as in our real IR- 
ACLEAN (see Section EH]). 

Figure [T7] shows the results of the Monte- Carlo sim- 
ulations. For all pairs with flux ratios > 40 and some 
pairs with separations < 2" 5, SExtractor cannot resolve 
them in the simulated J+Ks image. On such cases, IRA- 
CLEAN does not know there exist two objects, and thus 
we have no handle on them. These cases are not shown 
in Figure[l7l This is a limit of the WIRCam imaging, not 
IRACLEAN. On the other hand, when SExtractor is able 
to resolve the objects, IRACLEAN has very good perfor- 
mance for blended objects with separations greater than 
3'.'0 within the entire flux ratio range. For the fainter 
objects in pairs, the fluxes can be under- estimated by up 
to 50% with separations less than 3'.'0; the smaller the 
separation, the more severe the under-estimate is, and it 
is flux ratio dependent. For the bright objects in pairs, 
however, the fluxes can be over-estimated by up to 50% 
with separations less than 2'.'0 when flux ratios are less 
than 3.5. In these extreme cases, we are pushing both 
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Fig. 15. — IRACLEAN performance affected by the accuracy of the deconvolution PSF. The left panel is for IRAC 3.6 fim and the right 
panel is for IRAC 4.5 fim. The Y-axis is the ratio between the scatter of fluxes measured using different PSFs and the residual noise. The 
X-axis is the IRAC flux measured using the best PSF. The black solid line indicates a ratio of 1.0 for Y-axis The higher the flux rms to 
noise ratio is, the more serious the PSF effect is. The light grey data points indicate isolated objects defined as that they have no neighbors 
brighter than 30% of their fluxes within 6"0. The dark grey data points are the rest of objects. The ratio between the numbers of isolated 
and non-isolated objects is about 2:5. The IRACLEAN fluxes of the isolated and non-isolated objects can be affected by the accuracy of 
the deconvolution PSF because of different reasons. See text for details. 

IRACLEAN and SExtractor to their limits. 

To estimate how many sources in our catalog may suf- 
fer from the under-estimate issue of the fainter source 
in a pair, we used the results of the Monte-Carlo sim- 
ulations to determine the flux ratio (FR) vs. angular 
separation relation where the flux is under-estimated by 
20%, and the relation can be described by 



AS = 1.84 + 0.5 x log{FR), 



(2) 



where ^45* is the angular separation and FR is the flux 
ratio. In other words, if one pair with a flux ratio of 
FR has an angular separation less than AS, the flux 
of the fainter object in the pair is under-estimated by 
greater than 20% (but not by ^50% according to our 
simulations). We then used Equation [2] to estimate how 
many sources in our catalog are in the abovementioned 
situation. The numbers are 3347 (5.4%), 2840 (4.6%), 
1086 (1.7%), and 691 (1.1%) for 3.6/mi, 4.5/mi, 5.8/im, 
and 8.0/im, respectively. Fewer sources are effected in 
5.8/im and 8.0/im because of their much lower surface 
number densities due to the shallower detection limits 
(see Section |8~T1 for details). 

To estimate how many sources in our catalog may suf- 
fer from the over-estimate issue of the brighter object in 
a pair, we utilized Equation [2] with FR < 3.5 to repeat 
the counting but for the brighter sources. The numbers 
are 872 (1.4%), 707 (1.1%), 430 (0.7%), and 272 (0.4%) 
for 3.6/zm, 4.5/xm, 5.8/im, and 8.0/im, respectively. It is 
worth noting that the real number of objects which are 
affected by the under-estimate and over-estimate issues 
may be less than the abovementioned estimated num- 
ber because the FRs calculated using the fluxes of real 
objects in our catalog are already biased by the under- 
estimate and over-estimate issues so that the resulting 
AS' from Equation [5] are greater than the ideal one cal- 
culated using the unbiased fluxes. 

These results suggest that the performance of IRA- 
CLEAN is generally good and only a small number of 
objects in our catalog are affected by their neighbors, 



and the effect is not ^50%. This simple simulations, 
however, do not include different morphologies, number 
of blended objects, and S/N ratios. We therefore per- 
form another simulation using the mock IRAC images to 
take the abovementioned parameters into account, and 
show the results in Section 17.3.31 

7.3.3. Mock IRAC image simulations 

To evaluate the IRACLEAN performance affected by 
the effects of morphology, number of blended objects, 
and S/N ratio, We generated four mock IRAC images by 
convolving the J+K$ mosaic with the ideal IRAC PSFs 
for the four IRAC channels. We then ran IRACLEAN 
on these mock IRAC images with the segmentation map 
of the original J+K$ image. Assuming the SExtractor 
FLUX_AUTO measurement of each object in the orig- 
inal J+Ks image is the true flux amplitude 0, we can 
examine how well IRACLEAN performs under a compli- 
cated condition by checking if IRACLEAN can recover 
the SExtractor fluxes. The results are shown in Fig- 
ure [15] For 3.6 /im and 4.5 /im, the fluxes in general are 
not under-estimated according to their running medians. 
For 5.8 /im and 8.0 /im, however, the fluxes around 3.0 
arc under-estimated by 0.05 to 0.1 mag in general be- 
cause of their larger PSFs. We should note that for the 
real SIMPLE IRAC data, the depths in 5.8 /im and 8.0 
/im are much shallower than that in J, K$, 3.6 /im, and 
4.5 /im, so the surface number densities in 5.8 /im and 
8.0 /xm are much lower than that in the other bands. 
The under-estimated issue for 5.8 /im and 8.0 /im shown 
in our simulations should therefore be much milder in 

9 According to the SExtractor flags in the J+K$ catalog, 
about 70% of the total number of objects have significantly biased 
FLUX_AUTO because of bright and close neighbors, and/or were 
originally blended with another one. Therefore the FLUX_AUTO 
measurements for these cases many not be the true answer. How- 
ever, the mock IRAC images have even more severe blending issue 
as compared to the original J+Kg image. The assumption here 
therefore is still reasonable and feasible. 
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Fig. 16. — Performance of IRACLEAN. The SIMPLE IRAC images are shown in the left panels and their residual images are in the 
right panels. The region in each panel is 300" X 140". North is up and east is to the left. The brightness and contrast of each panel are 
identical. From top are the IRAC images of 3.6/^m (channel 1) to 8.0/^m (channel 4). All images are shown with inverted linear scales. 



the real SIMPLE data. For fluxes fainter than 2.0, they 
are over-estimated because of the selection effect in the 
fainter data. According to this plot, we conclude that 
the real SIMPLE IRACLEAN fluxes with greater than 
5cr detections would have biases <C0.1 mag. 

We also ran SExtractor in the single-image mode for 
the mock IRAC images and the results are shown in Fig- 
ure [19J According to Figure HH the under-estimate issue 
of SExtractor is more severe than that of IRACLEAN, 
and the scatters are at least a factor of 2 of that using IR- 
ACLEAN. The results suggest that IRACLEAN should 
perform better than SExtractor for the SIMPLE IRAC 
images. 

7.3.4. Comparison with the cross-convolution method 

In iWang et al. I (|2010l ). an alternative cross- 
convolution method (XCONV) was introduced for 
measuring ifs-IRAC color. The concept of this cross- 
convolution method is to match the PSF between K$ 
and each IRAC channel by convolving their PSFs to 
each other. Since this method does not require any 



IRACLEAN procedure but just simple photometric mea- 
surements on PSF-matched images, it should provide 
the most accurate ifs-IRAC colors for isolated bright 
objects. We thus compare our IRACLEAN results with 
XCONV colors. 

We first convolved each IRAC image with the 
WIRCam K s PSF, and convolved the WIRCam K s im- 
age with the IRAC PSFs we used in IRACLEAN. We 
take into account that there exist four different IRAC 
PSFs at 3.6 and 4.5 /mi (Section EU} ■ The FWHMs 
of the XCONVed PSFs are between 2" to 3". We then 
used the double-image mode of SExtractor to measure 
fluxes on these XCONVed images, by using the uncon- 
volved J+Ks image as the detection image. Because 
we focus on measuring Ks~lRAC colors here, we do not 
need to recover total fluxes since we have matched the 
PSFs in the Ks and IRAC images. We adopt an aper- 
ture size of 3" in diameter in all cases, to have a good 
balance between S /N and the inclusion of most fluxes for 
all kinds of source morp hologies. These are all similar to 
the XCONV method in IWang et aD (f2C)Toh . 
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Fig. 17. — Simple Monte-Carlo simulations of IRACLEAN for two blended objects. The simulation includes different angular separations 
and flux ratios between two blended objects. Four large panels indicate the simulation results for four IRAC channels. The FR value 
shown in each sub-panel is the flux ratio between two blended objects. The ratios between the output flux and the input flux of objects 
in the blended pairs are shown in the vertical axes. The open box indicates the brighter object in a pair while the filled box indicates the 
fainter one. Some missing filled boxes in each panel mean that SExtractor cannot successfully deblend the two objects in the WIRCam 
image; IRACLEAN is not able to handle this case since it relies on SExtractor for object detections. The IRACLEAN fluxes of the brighter 
objects in pairs therefore are the combined fluxes of both objects in pairs. 



The comparisons of the ifg-IRAC colors derived using 
the IRACLEAN and XCONV methods are shown in Fig- 
ure[5DJ The X5-IRAC colors derived using both methods 
are consistent with each other. The systematic offsets be- 
tween these two colors (XCONV-IRACLEAN) for 1,231 
objects with F Ks > 100/^Jy are -0.032, -0.011, -0.031, 
and -0.011 mag, for 3.6, 4.5, 5.8, and 8.0 ^m, respec- 
tively, which are within their statistical errors. We note 
that in the IRACLEAN case, the Ks fluxes are measured 
using the SExtractor AUTO aperture, which ca n miss up 
to 5% of total fluxes (|Bertin fc Arnouts |[l99(l . The IR- 
ACLEAN measurements of the IRAC fluxes, however, 
does not have such systematics. This may explain why 
the IRACLEAN colors are systematically redder than 
the XCONV colors. Nevertheless, the consistency be- 
tween the two colors is fairly good. We also checked the 



objects with color differences greater than 0.2 mag and 
Fks > 100/xJy. Most of them are in crowded areas with 
multiple bright sources so that their XCONV colors may 
be biased because of the large XCONV PSFs. Only a 
few of them are saturated/bright extended objects where 
their flux measurements are affected by the limitations 
of IRACLEAN, as we mentioned in Section ITXol 

7.3.5. Comparison with SExtractor 

Figure FJU shows a comparison of IRAC photometry 
using IRACLEAN against that using SExtractor with 
FLUX_AUTO in the single-image mode. The two pho- 
tometric methods are consistent on bright objects, but 
the scatter increases on fainter objects. Even on bright 
objects, as well as faint objects, there are small magni- 
tude offsets. The gray dashed lines in the figure indicate 
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Fig. 18. — IRACLEAN performance on the mock IRAC images. Four panels from left to right, top to bottom are for 3.6 /im, 4.5 (im, 5.8 
fim, and 8.0 fim. The Y-axis is the magnitude difference between the IRACLEAN flux in the mock IRAC image and the SExtractor flux in 
the original J+Kg image. The X-axis is the SExtractor flux in the original J+Kg image. The red line indicates the running median while 
the green lines indicate the upper and lower 68th percentiles of the distribution. The vertical grey line indicates the 5<r limit of IRACLEAN 
flux. 



a magnitude difference of -0.05, which well describe the 
differences in magnitudes derived using the two methods. 
This is consistent with th e 5% flux loss of SExtra ctor 
FLUX.AUTO discussed in lBertin fc Arnouts I (pM ). 

At 3.6 and 4.5 /im, the distributions of flux ratios are 
asymmetric about -0.05 mag. There are more objects 
whose SExtractor fluxes are larger than IRACLEAN 
fluxes. This is consistent with SExtractor fluxes being 
boosted by nearby objects. Such a trend is not apparent 
at 5.8 fim and is even reversed at 8.0 /Ltm. In these two 
bands, the sensitivity is lower and thus the mean sep- 
aration between detected objects is larger. This makes 
flux booting by nearby objects less an issue. In addition, 
the SExtractor auto aperture may miss a great portion 
of the flux of faint objects given the much broader PSFs 
at the longer IRAC wavelengths. 

Some bright sources have > 0.1 mag differences be- 
tween their IRACLEAN and SExtractor fluxes. Accord- 
ing to the flags provided by SExtractor, more than 70% 
of objects brighter than 50/jJy are blended with their 
neighbors. In most of these cases, the SExtractor fluxes 
are brighter than the IRACLEAN fluxes, consistent with 



their SExtractor fluxes being boosted by their neighbors. 
In a handful of cases, the large mag differences are a con- 
sequence of terminating IRACLEAN on saturated/bright 
extended objects as discussed in Section 17.2.61 We also 
verify that all the objects with large differences between 
their IRACLEAN and SExtractor fluxes are blended with 
their neighbors. According to Section 17X31 IRACLEAN 
can provide better flux estimates for faint sources in the 
IRAC images, and thus the majority of scatter in Fig- 
ure [2H would be due to the limitation of SExtractor. 

7.3.6. Comparison with the FIREWORKS Catalog 

We compared our IRA C fluxes with tho se from 
the FIREWORKS catalog (IWuvts et al. I l2008h . which 
provides the fluxes for the IRAC GOODS-S data. 
IWuvts et al. I (|2008l ) used the segmenta tion map gen- 
erated from the VLT ISAAC X s -band (jRetzlaff et al. I 
l2O~L0h image as a prior and convolved the -ftTs-band im- 
age with IRAC PSFs for each object. They then tried 
to remove the neighbors around a certain object in the 
IRAC images, by subtracting their flux-matched con- 
volved images. By assuming that all the neighbors are 
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Fig. 19. — Same as Figure [18] but using SExtractor. 

well-subtracted for the target, they then directly put 
an aperture to measure its flux. This is an alterna- 
tive method to minimize the blending issue for IRAC 
flux measurements, and it is very useful to compare the 
relative merits of IR ACLEAN and the method used by 
IWuvts et all (pOOl) . 

We show the comparison in Figure [2U In general the 
IRAC fluxes between the two catalogs are consistent with 
each other. The very small dispersions in Figure [H] es- 
pecially for 3.6 and 4.5 /xm, imply that these two entirely 
independent methods both work well in terms of resolv- 
ing the IRAC blending issue. There are, however, small 
differences between the two catalogs. At the bright end, 
objects brighter than 30 /xJy have ~ 5% higher fluxes 
in the FIREWORKS catalog. At the faint end, particu- 
larly at 5.8 and 8.0 fim, objects also have higher fluxes 
in the FIREWORKS catalog. The exact nature of these 
small differences are unclear to us. They are unlikely 
due to under-estimated IRACLEAN fluxes for blended 
objects, as we demonstrated in our Monte Carlo simula- 
tions (Section r7.3.2[) . W e suspect that t hey a re caused by 
the unmatched PSFs in IWuvts et al. I (|2008l) . since they 
did not deconvolve their Kg image with the Ks PSF be- 
fore the convolution with the IRAC PSFs, which may 
lead to over-subtracting the outer wings of the neighbor- 



ing objects and cause under-estimated local background. 
However, this is hard to verify with the dat a we have and 
witho ut knowing the exact procedure in IWuvts et al. I 
(|2008h . 

Notice that the scatter in the distributions at the faint 
end at 3.6 and 4.5 /im in Figure [55] are very small. This 
suggests that the reference image (i.e., the ISAAC Ks 
image) is not as deep as the IRAC 3.6 and 4.5 /im im- 
ages, so the fluxes of many fainter IRAC objects were not 
measured and included into the FIREWORKS catalog. 
Moreover, the faint IRAC objects that are not detected 
in the ISAAC Ks image were not subtracted before aper- 
ture photometry was applied to their brighter neighbors. 
Therefore, the fluxes of their brighter neighbors may be 
over-estimated. They may also contribute to the back- 
ground noise and make the FIREWORKS detection lim- 
its worse. 

7.3.7. Comparison with the SIMPLE Catalog 

IDamen et al. I (|2011D published a photometric catalog 
of the SIMPLE data. They used AUTO aperture in SEx- 
tractor to measure fluxes of all objects. If the AUTO 
aperture for an object is smaller than 4" in diameter, 
they replaced the AUTO flux with flux measured with 
a fixed 4" diameter aperture. If objects are blended, 
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Fig. 20. — Color differences between the IRACLEAN and XCONV methods. The red line indicates the running median and the green 
lines indicate the upper and lower 68th percentiles. The systematic offsets between these two colors for 1,231 objects with Fk s > 100 arc 
-0.032, -0.011, -0.031, and -0.011 for 3.6/xm, 4.5/rni, 5.8/zm, and 8.0/xm, respectively. The plot shows that the Kg - IRAC colors derived 
using both methods are consistent with each other. 



their 4"-aperture fluxes are also used. The y then applied 
apertu re corrections on all these objects. IDamen et al. I 
(|2011f ) excluded blended objects when they compared 
their source fluxes with the FIREWORKS catalog, claim- 
ing that these objects worsen the comparison. We thus 
first compare our fluxes with the SIMPLE catalog on iso- 
lated objects according to the blended flag in the SIM- 
PLE v3.0 catalog. The result is shown in Figure [2^1 and 
the comparison is fairly good. We emphasize, however, 
that more than 70% of the matched objects are marked 
as blended sources in the SIMPLE catalog. The result 
shown in Figure [23] is therefore not a representative com- 
parison between the two catalogs. 

We next compare our IRACLEAN results with the 
SIMPLE v3.0 catalog on blended objects. The results, 
shown in Figure 1241 are quite striking. There are two 
distinct sequences in each IRAC band. The upper se- 
quences contain approximately 66%, 68%, 54% and 51% 



of all blended objects at 3.6, 4.5, 5.8, and 8.0 /im, re- 
spectively. The SIMPLE fluxes of the upper sequences 
are only slightly higher compared with our IRACLEAN 
fluxes. The reason for the differences is that the SIM- 
PLE fluxes of blended sources would be over-estimated 
because conventional aperture photometry method was 
used, while the IRACLEAN method is designed to esti- 
mate relatively unbiased fluxes for blended sources. On 
the other hand, the SIMPLE fluxes of the lower sequences 
are between 30% to 60% lower than the IRACLEAN 
fluxes. This very large offset, and the fact that there 
are two distinct sequences, cannot be explained by any 
systematic effects that we are aware of. We looked at 
the objects in the lower sequence objects in the IRAC 
images, but did not find anything different about these 
objects compared with those in the upper sequences. We 
remind the reader at thi s point that both IR ACLEAN 
and the method used by IDamen et al. I (|2011D to derive 
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Fig. 21. — The comparison of the IRAC photometry between the IRACLEAN and FLUX_AUTO methods. X-axis is the flux for each 
IRAC channel and Y-axis is the magnitude difference (or flux ratio) between IRACLEAN and FLUX_AUTO methods. The black line 
indicates offset while the gray dashed line indicates an offset of -0.05 for Y-axis. The red line indicates the running median and the green 
lines indicate the upper and l ower 68th percentile s . The offsets between the results derived using the two methods are consistent with the 
missing flux ratio provided by Bcrtin & Arnouts (1996). 



fluxes of indiviaul objects are based on the same SIMPLE 
IRAC images, yet we do not see such second sequences 
when we compare our results with XCONV, SExtractor, 
and FIREWORKS. The exhaustive comparisons that we 
have made with all relevant existing catalogs suggest that 
the lower sequences are previously unknown systematic 
effects in the SIMPLE v3.0 catalog. We are therefore 
concerned about the conclusions reached in studies that 
have been made based on the SIMPLE v3.0 catalog. 

8. COMBINED NIR AND IRAC CATALOG 

We combined the TENIS WIRCam photometry and 
the SIMPLE IRACLEAN photometry. We also removed 
objects with fluxes less than 3er in both the J and Kg 
bands. The final TENIS WIRCam and IRAC catalog 
includes 62,326 objects. We did not apply Galactic ex- 
tinction correction to this catalog. We also include the 
XCONV colors since they are good color references for 



isolated bright sources. Table [2] shows the statistics of 
this catalog which has the following format: 

(1) : Object ID. 

(2) - (3): R.A. and Decl. 

(4) — (15): Flux and flux error measurements (in /iJy) 
for J, Kg, IRAC 3.6/im, 4.5/im, 5.8/xm, and 8.0/irn. 

(16) - (19): XCONV K s - IRAC colors (in AB mag- 
nitude) . 

8.1. Depths 

The sensitivity of our catalog are by far the highest 
in all the wavebands from J to 8.0 fim among exist- 
ing observations in the ECDFS and GOODS-S. In the 
ECDFS region, the median 5a limits among all detected 
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Fig. 22. — IRAC photometry comparison between our catak 
indicates the running median while the green lines indicate the 
fluxes in the two catalogs with an rms of flux ratio dispersions o) 

TABLE 2 
Number of Objects and Detection 
Limits in the Multiband Catalog 



Waveband Detections Depth (nJy) 



J 53722 0.115 

K s 57492 0.199 

3.6,um 53801 0.083 

4.5A*m 49667 0.105 

5.8tim 28284 0.541 

8.0/xm 22418 0.669 



Note. — Detections arc the number 
of objects with S/N > 3. Depths are the 
medium values of lcj flux errors for objects 
with S/N > 3. 

objects in our catalog are 24.50, 23.91, 24.85, 24.60, 
22.82, and 22.59 mag for J, K s , 3.6, 4.5, 5.8, and 
8.0 /xm, respectively. These J and Ks limiting magni- 
tudes are much deeper than those of the MUSYC survey 



ig (IRACLEAN method) and the FIREWORKS catalog. The red line 
upper and lower 68th percentiles. The majority of objects have similar 
'. less than 8%. 

(jCardamone et al. 2010|). As shown in Sectional our K$ 
depth in the GOODS-S region is also ~ 0.5 mag deeper 
than t he VLT ISAAC data published in iRetzlaff et al~l 
(2010). The same is also true for our J depth, as com- 
pared to the VLT ISAAC J depth in iGrazian et al. I 
(200i). 

We found that our IRAC limitin g magnitudes a r e muc h 
deeper than those provided in ID amen et al. I (|2011D . 
which are 23.8, 23.6, 21.9, and 21.7 mag at 3.6, 4.5, 
5.8, and 8.0 /zm, repectively. Since our work and that 
of Damen et al. are based on the same IRAC dataset, 
it is thus important to examine whether our sensitivities 
are reasonable. We first compare the residual noises in 
our IRACLEAN images with those provided by Spitzer 
Sensitivity Performance Estimation Tool (SENS-PET). 
The detection limit estimated by SENS-PET is derived 
using a 10-pixel radius aperture (~ 24" in diameter) for 
point-like sources. We therefore measured the rms of 
the background noise in our IRACLEAN residual images 
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Fig. 23. — Photometry comparison between the IRACLEAN flux and the SIMPLE catalog. The red line indicat es the running media n 
while the green lines indicate the upper and lower 68th percentiles. Only isolated objects are shown as suggested by Damcn e t al. I l|2011h . 
According to this figure, the fluxes of isolated objects are consistent between both catalogs very well. 



with a 24" diameter aperture. They are 0.28, 0.40, 1.74, 
and 1.74 /j,3y for 3.6, 4.5, 5.8, and 8.0^m, respectively, 
and where the average integration time in the SIMPLE 
IRAC images is - 1.5 hr. According to the SENS-PET, 
1 hr of integration will provide la sensitivities of 0.191, 
0.277, 1.56, and 1.67 fjjy for 3.6, 4.5, 5.8, and 8.0/xm, 
respectively. Our sensitivities are worse, which is ex- 
pected since the sensitivities quoted by the SENS-PET 
correspond to the ideal cases. Next, we convert the 
SENS-PET sensitivities to our IRACLEAN 2'.' 7 x 2'.' 7 
aperture, by assuming that the photometric error scales 
with square-root of the aperture area. After taking into 
account the aperture corrections ( Section I7.2.4[) . we ob- 
tained 5 a limiting magnitudes of 25.51, 25.12, 23.12, 
and 22.93 for 3.6, 4.5, 5.8, and 8.0 ^im, respectively. Our 
actual limiting magnitudes are again shallower. 

The above comparisons show that our sensitivities are 
better than what had been previously achieved on the 
SIMPLE IRAC data, but do not exceed what can be 



achieved in ideal cases. The differences in the detec- 
tion li mits between our catalog and that of lDamen et al. I 
(| 2 1 If) are most likely caused by confusion effects in the 
IRAC images. With our J+K$ prior image that is nearly 
as deep as the IRAC images, IRACLEAN is much less 
affected by faint undetected sources as well as the PSF 
wings of nearby bright objects. This allows us to esti- 
mate fluxes for faint sources more accurately, and thus 
pushing the detection limits closer to the ideal values. 

8.2. Spurious sources 

As we mentioned in the beginning of this section, we 
cleaned the TENIS multi-wavelength catalog by remov- 
ing objects with low S/N. Although this step should have 
eliminated most spurious objects, the catalog could still 
be significantly contaminated. We therefore investigated 
the spurious fraction in our catalog. First, we inverted 
all the images (i.e., making negative images) including 
the J, Ks, J+Ks, and the four IRAC images. We 
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Fig. 24. — Same as Figure l23l but for blended objects. 

then repeated identical steps used to generate the multi- 
wavelength catalog on these inverted images. The final 
catalog for the inverted images contains 21578 objects, 
which suggests that the spurious fraction of the TENIS 
multi- wavelength catalog is unreasonably high (~ 35%). 
We checked the inverted J+K$ image and the "sources" , 
and we found that most of the "sources" correspond to 
negative holes produced by the crosstalk removal proce- 
dure of the WIRCam reduction pipeline. (The crosstalk 
removal does not generate posit i ve fea tures.) This is sim- 
ilar to the case in lWang et al. I (j2010). If we just use the 
relatively crosstalk-free regions to calculate the spurious 
fraction, the value dramatically decreases to 6%. We 
note that the value of 6% is an upper limit since fainter 
objects still can produce low-level crosstalk. 

A second test on the spurious fraction is objects de- 
tected by both WIRCam and IRAC. According to Ta- 
ble [21 there are 53801 objects detected at 3.6/Ltm with 
S/N > 3, suggesting a spurious fraction of 14%. How- 
ever, this is also an upper limit. The J and 3.6/im color- 
magnitude diagram in Figure [25] shows that blue galaxies 



with J - 3.6/im < -0.4 start disappearing at J > 25. This 
indicates that many IRAC undetected but J detected ob- 
jects are real. Therefore, the spurious fraction must be 
much less than 14%. Based on the above two tests, we 
conclude that the spurious fraction of the TENIS catalog 
is less than 6%. 

We do not attempt to quantify the completeness of our 
catalog, which is nontrivial given the complex nature of 
IRACLEAN and the fact that we detect sources from a 
J+Ks image. Readers interested in using our catalog 
and wishing to know the completeness at a given flux 
level should run their own source extraction and selec- 
tion, quantify the completeness, and then extract pho- 
tometry from our multiwavelength catalog. 

9. SUMMARY 

We present an ultra-deep J and Ks dataset covering 
a 30' x 30' area in the ECDFS, as part of our Taiwan 
ECDFS Near-Infrared Survey (TENIS). The median 5a 
limiting magnitudes for all objects reach 24.5 and 23.9 
mag (AB) for J and Ks, respectively. In the inner 400 
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IRAC fluxes for most objects. We also compared IR- 
ACLEAN fluxes with fluxes directly measured by SEx- 
tractor, and with the FIRWORKS and SIMPLE cata- 
logs. We found that IRACLEAN results are superior in 
many cases, and our IRACLEAN results provide by far 
the deepest IRAC catalog in the ECDFS region. This 
J+Ks detected catalog consists of flux measurements 
for J, K s , IRAC 3.6, 4.5, 5.8, and 8.0 fim, and XCONV 
K S -IRAC colors for all four IRAC bands. We publicly 
release the data products of this work, including the J 
and Ks images, and the J+Ks selected multiwavelength 
catalog. 



Fig. 25. — J and 3.6/^m color-magnitude diagram. The gray 
solid-line indicates the 5<x limiting magnitude for point sources in 
J, and the gray dashed-line indicates the 5<r limit of J-3.6/im color 
for point sources. The blue galaxies with J — 3.6/im < —0.4 start 
disappearing at J > 25 (the bottom-right corner), which implies 
that many objects in the TENIS catalog with S/N < 3 at 3.6/im 
are real J-detectcd objects. 



arcmin 2 region of the images where the sensitivity is more 
uniform, objects as faint as 25.6 and 25.0 mag are de- 
tected at 5 a. In our final catalog, we detect objects in 
a J+Ks image in order to achieve higher completeness. 
We also developed a novel deconvolution technique (IR- 
ACLEAN) to accurately estimate the IRAC fluxes for 
all the J+Ks detected objects in the ECDFS, using our 
J+K image as a prior. With simple Monte- Carlo simu- 
lations and comparison against the XCONV technique, 
we showed that IRACLEAN is able to correctly recover 
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